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CN ■ Abstract 
> , 

! We study the linear perturbations of multi-field inflationary mod- 

els governed by a Lagrangian which is a general function of the scalar 
fields and of a global kinetic term combining their spacetime gradients 
with an arbitrary field space metric. Our analysis includes ^-inflation, 
Dirac-Born-Infeld inflation and its multi-field extensions which have 
been recently studied. For this general class of models, we calculate 
the action to second order in the linear perturbations. We decompose 
the perturbations into an adiabatic mode, parallel to the background 
^ ' trajectory, and entropy modes. We show that all the entropy modes 

propagate with the speed of light whereas the adiabatic mode prop- 
agates with an effective speed of sound. We also identify the specific 
combination of entropy modes which sources the curvature perturba- 
tion on large scales. We then study in some detail the case of two 
scalar fields: we write explicitly the equations of motion for the adia- 
batic and entropy modes in a compact form and discuss their quantum 
fluctuations and primordial power spectra. 
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1 Introduction 



Inflation has now become a standard paradigm for describing tlie pliysics of 
the very early universe, but the nature of the field(s) responsible for inflation 
remains an open question. The hope is that future cosmological observations, 
in particular those of the CMB, will be able to rule out large classes of models 
and give some hints on the underlying physics. 

In the last few years, intensive effort has been devoted to trying to con- 
nect string theory and inflation (for recent reviews, see e.g. [I]- [5]). For 
simplicity, most studies of string inflation have considered a single effective 
scalar field. However, in the low-energy limit of string theory, many scalar 
fields are present and several of them could thus play a dynamical role during 
inflation. This would affect the generation of primordial perturbations. For 
instance, whereas in single-field inflation, the curvature perturbation is con- 
served on large scaled], the curvature perturbation in multi-field inflation can 
generically be modified on large scales, because it is sourced by entropy (or 
isocurvature) perturbations. This feature was first pointed out in [S] in the 
context of Jordan-Brans-Dicke type gravity where the gravitational sector 
contains a scalar field. This effect has also been illustrated recently P, [10] in 
the context of specific inflationary models based on string theory construc- 
tions, and shows that the restriction to an effectively single-field scenario, 
despite its appealing simplicity, might be misleading as the final curvature 
perturbation, which will be eventually observed, can sometimes originate 
mainly from the entropy modes. 

In multi-field inflation, there is also the possibility, depending on the re- 
heating scenario, to produce, after inflation, both adiabatic and isocurvature 
perturbations, which can be correlated, as first pointed out in [TT] . The CMB 
measurements have shown that the primordial perturbations are mainly adi- 
abatic but a small amount of isocurvature modes is still allowed by the data 
(see e.g. [I2]-[T5] for the most recent analyses). 

In this context, a lot of works have been devoted to multi- field inflation. 
Although many of these studies usually assume standard kinetic terms for 
all the scalar fields involved, several consider multi-field inflation with non- 
standard kinetic terms described by a non-trivial metric Gij{(j)^) in field 
space, so that the kinetic part of the Lagrangian is of the form 

Cmu = -\Gijd,(t>'d^ct>'. 

More sophisticated types of non-canonical kinetic terms can also be envis- 
aged. This is in particular the case for an interesting scenario based on string 

^see however [S] and [7] for exceptions. 
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theory, in which the inflaton is a scalar field characterizing the position of 
a probe brane moving in a warped background geometry. This scenario has 
been named DBI (Dirac-Born-Infeld) infiation [T6l [T7] , as the infiaton field 
is governed by a Dirac-Born-Infeld action. The latter can be seen as a par- 
ticular case of the more general framework of infiation [181 CH] 5 where the 
action is an arbitrary function of the infiaton and of the square of its space- 
time gradient. Two recent works [201 EI] have studied the perturbations in 
effectively multi-field DBI infiation, where the extra fields correspond to the 
angular degrees of freedom of the moving brane |22j-|25j. 

The purpose of the present work is to analyse a very large class of multi- 
field models, which can be described by an action of the form 



where P is an arbitrary function of scalar fields and of the kinetic term 



where G/j = Gij{(j)^) is an arbitrary metric on the A^-dimensional field 
space. This can be seen as a generalization of the Lagrangian of /c-infiation 
|18j to the case of several scalar fields. This also includes DBI infiation and 
its multi- field extensions studied in [20l [21]. 

In this work, we expand the above action up to second order in the 
linear perturbations. The second order action can be used to derive the 
classical equations of motion for the perturbations. It is also the starting 
point to calculate the spectra of the primordial perturbations generated from 
the vacuum quantum fiuctuations of the scalar fields during infiation. 

In order to analyse the equations of motion, we decompose the linear 
perturbations into the (instantaneous) adiabatic perturbation, i.e. the per- 
turbation parallel to the background trajectory in field space, and the (in- 
stantaneous) entropy perturbations, which are orthogonal (with respect to 
the field space metric) to the trajectory. We find, quite generically, that the 
equation of motion for the adiabatic perturbation is a wave equation which 
depends on an effective sound speed Cg while all the entropy perturbations 
obey a wave equation involving the speed of light. This shows that the prop- 
erty that adiabatic and isocurvature perturbations propagate with different 
speeds, pointed out in [20] for the particular case of two-field DBI infiation, is 
valid for the large class of multi-field models studied here, whatever the spe- 
cific dependence of the action on the kinetic term, and whatever the number 
of scalar fields involved. For more than two scalar fields, the entropy sec- 
tor contains several degrees of freedom but they all satisfy a light-like wave 
equation. 




(1) 



(2) 
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We also analyse the evolution of the curvature perturbation on large 
scales, and show that it is sourced by a specific combination of entropy pertur- 
bations, which can be decomposed into a term corresponding to the bending 
of the background trajectory in field space, as in multi-field infiation with 
canonical kinetic terms, and an additional term, which is present only for 
non-trivial functions P. We finally specialize our analysis to the case of two 
fields and compare our results with previous works. 

The plan of this paper is the following. In the next section, the back- 
ground equations of motion are given. Section 3 is devoted to the derivation 
of the second-order action for the linear perturbations. In the subsequent sec- 
tion, we analyse the (comoving) curvature perturbation and its large-scale 
evolution. In Section 5, we restrict our analysis to two scalar fields and 
study in detail the equations of motion for the adiabatic and isocurvature 
components. Using the second-order action, we then discuss, in Section 6, 
the quantization of the perturbations in the simple case where the adiabatic 
and isocurvature modes are decoupled. We conclude in the final section and 
present in an appendix the application of our formalism to the DBI case. 

2 Background 

As discussed in the introduction, our starting point is the action 

S = |d^xv^(ii? + P(X,0O) (3) 

with 

X = -^GjjV^<P'V>'<P'. (4) 

where we have set SttG = 1 for simplicity. Throughout this paper, we will 
use the implicit summation rule on the field indices I, J, . . .. 

The energy-momentum tensor, derived from ([H]), is of the form 

T^"" = Pg^'" + PxGijVy^Vy-^ , (5) 

where Px is the partial derivative of P with respect to X. The equations of 
motion for the scalar fields, which can be seen as generalized Klein-Gordon 
equations, are obtained from the variation of the action with respect to 0^. 
One finds 

"^.{PxGjjV^q^') - ^Px(V^0^)(V>^)9,GxL + Pi = 0. (6) 
where P/ denotes the partial derivative of P with respect to (p^ . 
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In a spatially flat FLRW (Friedmann-Lemaitre-Robertson- Walker) space- 
time, with metric 

ds^ = -dt^ + a?{t)dx^, (7) 

the scalar fields are homogeneous and the energy-momentum tensor reduces 
to that of a perfect fiuid with energy density 

p = 2XPx-P, (8) 

and pressure P. 

The evolution of the scale factor a{t) is governed by the Friedmann equa- 
tions, which can be written in the form 

if2 = i(2XPx-P), (9) 

and 

H = -XP,x. (10) 
The equations of motion ([6]) for the scalar fields reduce to 



+ T'j^^U'' +{3H + ^U'- -^G'^Pj = , (11) 



where the Fj^ denote the Christoffel symbols associated with the field space 
metric Gjj. By noting that the first two terms in the above equation simply 
correspond to the components of the acceleration in curved coordinates (here 
in field space), which we can write as 

A0'^0' + r^j,,0^0^, (12) 

equation (ITT]) can be rewritten as 

+ + ) 0' - -^^''^^J = ' (13) 
or, in an even more compact form, as 

a-^Vt{a'Px<f>i)=Pi, (14) 

where we have used the field metric Gjj to lower the field index J, so that 
<pj = Gjj(p'^; Dt acts as an ordinary time derivative on field space scalars (i.e. 
quantities without field space indices) and VtGjj = 0. 
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The above expressions are compact but second derivatives of the scalar 
fields are hidden in the term Px- It is sometimes more useful to express the 
equations of motion in the form 




P,K^A^' -^G''P,j = Q. (15) 



where we have introduced, as in [T9] . 




P,x 



(16) 



which can be interpreted as the square of an effective sound speed, as we will 
see later. We will also use the dimensionless parameter 



3 Dynamics of the linear perturbations 

We now study the linear perturbations about the background solution dis- 
cussed in the previous section. A priori, one must consider the perturba- 
tions of the scalar fields as well as the metric perturbations. The scalar 
metric perturbations are coupled to the scalar field perturbations via the 
scalar field equations of motion and Einstein's equations. One could write 
down directly the linearized version of these equations. However, since we 
will be interested by the quantum fiuctuations of the perturbations during 
the infiationary phase, it is useful to compute the action at second order in 
the perturbations. From this action, one can easily derive the equations of 
motion for the linear perturbations. One can also determine the normaliza- 
tion of the vacuum quantum fiuctuations, and therefore the amplitude of the 
primordial cosmological perturbations. 

In order to quantize the system, it is useful to write the second-order ac- 
tion as an expression depending only on the true physical degrees of freedom. 
This can be done by using the constraints to simplify the second-order action 
but this is a tedious procedure [2S]- A quicker method uses the constraints, 
within the Hamiltonian formalism, as Hamilton- Jacobi equations in order to 
identify automatically the physical degrees of freedom [27] (see also [2H])- An 
even simpler procedure, introduced in [29], consists in writing the action in 
the ADM form (20] and in solving explicitly the constraints for the lapse and 
the shift. 

In the case of a single scalar field, two gauge choices are natural: the first 
is the gauge where the scalar field is spatially uniform on constant time slices; 



c. 



■s 



(17) 



s = 



6 



the second is the spatially fiat gauge, where the spatial part of the metric is 
unperturbed. The physical perturbation is described only by the metric in the 
first case, and only by the scalar field perturbation in the second case. When 
dealing with several scalar fields, the first possibility is no longer possible 
and the only natural choice is to go in the spatially fiat gauge, where the 
scalar field perturbations correspond to the physical degrees of freedom. The 
calculation which follows can be seen as a generalization (and a unification) 
of the second-order actions computed in [HI [32] with the same procedure. 



3.1 ADM formalism and constraints 

In the ADM formalism [33], the metric is written in the form 

ds"^ = -N^dt^ + hij {dx' + N'dt) {dx^ + N^dt) (18) 

where is the lapse function and A^* the shift vector. The action (I3l) then 
reads 

S = ]^j dtd^'xVhN (i?(=^) +2P) + ^J dtd^x^{EijE'' - E^) , (19) 

where h = det{hij) and R^^^ is the Ricci curvature calculated with hij. The 
symmetric tensor Eij, defined by 

Ei, = - A^(.b) (20) 

(the symbol | denotes the spatial covariant derivative associated with the 
spatial metric hij), is proportional to the extrinsic curvature of the spatial 
slices. The function P in fll9p depends on the kinetic term X, which can be 
written as 

^ = - ^h^'d4'd,ct>' (21) 

with 

= <j)^ - N^dj<f)^ . (22) 
The variation of the action with respect to A^ yields the energy constraint, 

2P - -^{E,,E'^ -E^ + 2P,xGijv'v') = , (23) 

while the variation of the action with respect to the shift A, gives the mo- 
mentum constraint, 

l^{Ei-E6i)y = ^Guv'd,<l>'. (24) 
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In order to study the linear perturbations about the FLRW background, 
we now restrict ourselves to the flat gauge, which corresponds to the choice 



hij = a\t)5ij. (25) 

The scalar fields on the corresponding flat hypersurfaces can be decomposed 
into 

0' = 0o + Q', (26) 

where the c^g ^^e spatially homogeneous background values of the flelds 
and the represent the linear perturbations. In the following, we will usu- 
ally drop the subscript '0' on 0q and simply identify cf)^ with the homogeneous 
background flelds. 

We also can also write the (scalarly) perturbed lapse and shift as 

iV = 1 + «, iVi = A* , (27) 

where the linear perturbations a and (5 are determined in terms of the scalar 
field perturbations by solving the linearized constraints. At first-order, 
the momentum constraint implies 



a 



2H 

while the energy constraint yields 



W , (2^ 



(29) 

where we have extended the notation P^, introduced in (fT2!) . to Q^, so that 

VtQ' = Q' + T'jK^'Q''. (30) 



3.2 Second-order action 

We now expand, up to quadratic order, the action in terms of the linear 
perturbations. This action can be written solely in terms of the physical 
degrees of freedom by substituting the expressions fpHjl and fpUl) for a and 
(3. In fact, it turns out that (3 disappears of the second order action, after an 
integration by parts. The calculation is straightforward, although somewhat 
lengthy. It is also convenient to regroup the terms involving derivatives of 

^In the ADM formalism, it is sufficient to use the perturbed lapse and shift up to first 
order as their second-order parts cancel out in the action. 
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the field space metric Gjj into covariant expressions. For example, we use 
the identity 



dt d^x 



dtd^xa^ 



P^xVtQ'GijVtQ' + {ViVjP + Px7^/xLJ0■^0^) Q'Q'] , (31) 



where we have discarded total derivatives. denotes the covariant deriva- 
tive associated to Gu (we thus have VjVjP = Pjj — T^jP^k) and TZijkl is 
the Riemann tensor for Gjj. 

The second-order action can be finally written in the rather simple form 



(2) 



dt d X a 



PxGrj + p. 



xx(Pi(PJ 



\x 



Gijd.Q'd'Q' 



-MjjQ'Q' + 2PxjhQ'VtQ[ 
with the effective (squared) mass matrix 

Mij = 



(32) 



-VjVjP - PxTZiKLJ^ 
1 . : : 



H 



xj9i 



P 



XI(PJ) 



-1(1 



2/^2 a 
The last term can be expanded, which yields 

XP^x 













4>i4>j 



.(33) 



Mij = 



p2 



-ViVjP - P^x'RiklA''^'^ + 



H 



P 



XJ<Pl 



+ P. 



xi(Pj 



XP\ sP\ 3P^x 



1 \ 



2H 



t<Pj 



2 



P 



X 



H 



1 + ^1 [P,xxX + P 



.XKS 



(34) 



The second-order action fl32l) is one of the main results of the present work 
and generalizes more restrictive cases considered in the literature [SU [2H1 E2] • 
Let us now discuss some of its features. First, one can notice that the non- 
fiat nature of the field space metric Gjj manifests itself in the replacement of 
the ordinary time derivatives by covariant time derivatives Vt and of Pjj by 
J)jT>jP, as well as the presence of the Riemann tensor TZijkl in the matrix 

Mjj [MIES]. 

A second interesting consequence of fl32l) is how the non-trivial depen- 
dence of the initial action on the kinetic term X affects the dynamics of 
the perturbations. One sees that the term involving the spatial gradients 
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is simply rescaled by a factor Px with respect to the standard case (where 
P = X — V so that Px = 1), which multiphes the metric Gjj. The term 
quadratic in the time derivatives is changed in a subtler way, since instead 
of P^xGij, one finds 

PxGjj = Px (gu + ^<i^i<Pj^ ■ (35) 

This shows that the background velocity, of components 0^, represents a 
special direction in field space as far as the dynamics of the perturbations 
is concerned. This direction corresponds to the (instantaneous) adiabatic 
direction, which has been introduced in [STj for multi-field inflation with 
standard kinetic terms. Introducing the 'adiabatic' unit vector e(, defined as 



ei = 4^, (36) 



and using (fTB|) . one finds 



Gjj^\e]e]+(Gjj-e]e]j), (37) 



where the term between parentheses represents, in field space, the projection 
orthogonal to the adiabatic direction. This decomposition clearly shows that 
the adiabatic part of the perturbations, i.e. along the field velocity, obeys 
a wave equation where the propagation speed is determined by the sound 
speed Cs, while the entropy perturbations, i.e. orthogonal to the scalar field 
velocity, propagate with the speed of light. This property, which was pointed 
out in [2ni EI] for the specific case of two-field DBI inflation, turns out to be 
generic for any system governed by an action of the form ([3]). 

The decomposition into adiabatic and isocurvature modes can be made 
more explicit by introducing an orthonormal basis in field space, {e^} [n = 
1, . . . , A^), where the first vector is the unit adiabatic vector introduced in 
flHUj) . The — 1 remaining vectors thus span the entropy subspace, which 
is orthogonal to the adiabatic direction. The decomposition of the perturba- 
tions on this new basis reads 

Q' = Q"e^, (38) 

with an implicit summation on the index n. In order to replace the time 
derivatives of the in terms of the time derivatives of the new components 
Q", one needs to take into account the time derivative of the basis vectors. 
It is in particular useful to define the quantities 

Zmn = CmlVtel^, (39) 
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which satisfy the antisymmetry property = — ^nm as a consequence of 
Vt{e^eni) = 0. The second order action then reads 



^2 = ^Jdtd^xa^ 



p 



a2 



'(40) 



'^s n>l 



where Mnm = ^li^u^m are the components of the interaction matrix in our 
new basis. Similarly, P^xn = ^i,P,xi and P„ = e^P/. 

So far, only the first vector of the basis, ei, has been specified. In the 
following, we will also specify the second element of the basis, 63, as the 
unit vector pointing along the projection on the entropy subspace of the 
field acceleration Vt(f)^ . The (N — 2) remaining vectors of the basis are left 
arbitrary in the present work. Note that [M] adopted a specific choice for all 
the vectors of the basis by considering the successive time derivatives of the 
background scalar fields. 



4 Curvature perturbations 

In the previous section, we have obtained the second order action in terms 
of the physical degrees of freedom . This completely determines the full 
dynamics of the linear perturbations. It is then useful to relate the quantities 

to other perturbed quantities which can be of interest, in particular the 
gauge-invariant quantities describing the metric perturbations. 

Instead of the metric written in the ADM form, we will work in this section 
with the (scalarly) perturbed FLRW metric written in the usual form (see 
[381126] for detailed reviews on the theory of linear cosmological perturbations 
and e.g. for a pedagogical introduction), 

ds"^ = -(1 + 2A)df + 2adiBdx'dt + [(1 - 2ij)Sij + 2dijE] dx'dx^ . (41) 

These metric perturbations can be combined to give the familiar gauge- 
invariant Bardeen potentials, defined by 



d 



a\E-B/a) , (42) 



$ = A 

dt 

* = ^ + a^H{E - B/a) . (43) 

For the matter, the linear perturbations of the energy density, pressure 
and momentum follow from the linearized energy-momentum, according to 
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the expressions 



Sp^-ST^, 6P = ^6Ti, dM = 5Tl (44) 



while, in our case, the anisotropic stress vanishes. The hnear combination 

= 5p- 3HSq (45) 

is gauge-invariant and defines the so-called comoving energy density pertur- 
bation (i.e. coincides with the energy density perturbation in the comoving 
gauge, characterized by 6q = 0). 

Other useful gauge-invariant quantities are obtained by combining the 
metric and matter perturbations. One can define the curvature perturbation 
on uniform-density hypersurfaces 

-C = ^P + -Sp, (46) 
P 

and the so-called comoving curvature perturbation 

n^^-^-Sq. (47) 

In the following, we will work in Fourier space for all these linear perturba- 
tions. 

The gauge-invariant quantities defined above can be related to the physi- 
cal degrees of freedom , which we have chosen to describe our system. The 
relations follow from the linearized Einstein's equations SG^j^iy = ST^^, more 
precisely from the constraints. In particular, the energy constraint reads 

3H{ij + HA) + —[ij + H{a^E - aB)] = --5p, (48) 

while the momentum constraint yields 

iP + HA = -^5q. (49) 

The combination of these two constraints gives a gauge-invariant relativistic 
generalisation of the Poisson equation, 

^vl/ = -l5p^. (50) 

Note that the Q^'s correspond to the gauge-invariant combinations = 
6(j)^ + {(f)^ /H) ip where the (50^'s are the field perturbations in any gauge. In 
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order to write the various quantities introduced above in terms of the , it 
is useful to notice that the ADM metric is equivalent to a metric of the form 
(l4T]l with the identification A = a, aB = f], ip = and ^ = 0. Using (l28ll . 
one finds that the perturbed energy density and pressure, in the flat gauge, 
take the form 



'flat 



Fx 
Hcl 



ljVtiHQ') + {2XPxi-Pi)Q' 



and 



SB 



fiat 



H 



bjVt{HQ') + PjQ'. 



(51) 



(52) 



Moreover, substituting (1281) in ( H9l) yields the expression of Sq, in the flat 
gauge, in terms of the . One thus gets 



and 



P 



X 



'^jVtQ' + 2XP^xiQ' - P,iQ' 



Px XP 



X 



H 



(53) 



?>hA(PiQ'. (54) 



Taking the time derivative of the expression (l53i) for TZ and combining the 
result with fl54p and (l5Up . one obtains an important result for the time evo- 
lution of IZ: 



n 



H a2 



H 



2XP 



X 



;i + cl)Pj - 2Xc^PiJ 



(55) 



where, for any covector Aj in field space, Af = Aj — (ef-ylx)ei/ represents 
its projection on the isocurvature subspace, i.e. orthogonal to the adiabatic 
direction. 

All the above results can be rewritten in the basis {e^}, which distin- 
guishes the adiabatic and isocurvature perturbations. One can express the 
comoving curvature perturbation as 



n 



H 



'2X 



(56) 



which illustrates that TZ characterizes the purely adiabatic part of the per- 
turbations. The comoving energy density perturbation now reads 



5p„ 



PxV2X 



H 2X) 



-2X P,XnQ' 
\n>l 



P,2Q'. 

(57) 
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Finally, the time evolution of TZ can be rewritten as 



• Hc'',k\ H 

TZ = ^ — \& H 

H 2XP 



X L n>l 

When one can ignore the first term on the right hand side, on sufficiently 
large scales, one recovers the familiar result that the curvature perturbation 
is sourced by entropy perturbations only. However, in our case, this entropy 
source term contains two contributions. 

The first one, proportional to P2, is a generalization of the term obtained 
in [37] . Indeed, for P = X — V{(j)^) and a fiat field space metric, the above 
relation with two fields reduces to 

^ = 1^ ^^g^ = ^-^"^+—668 , {P = X-V,Gu = M (59) 

where we have introduced, in the second equality, a = V2X, 5s = and 
9 = —V2/0', corresponding to the notation of [HT] . 

This result has been extended in [36] to the multi-field situation, still 
for P = X — V but with a general metric on field space. In this case, the 
curvature perturbation on large scales is still sourced by a term related to the 
perturbation in the direction 63, i.e. along the entropic projection of the field 
acceleration Vt(f)^ . Therefore, on large scales, the curvature perturbation is 
not conserved if the trajectory in field space is not geodesic. By a geodesic 
trajectory, we mean here that T>t4>^ is proportional to <p^ . Equivalently, in 
view of the background equations for the scalar fields ( |T3l) . the curvature 
perturbation is not conserved on large scales if the gradient of P along 63 
does not vanish. 

In the more general situation considered here, one sees that the mixed 
derivatives P^xi along the entropy directions give an additional contribution 
to the non- conservation of the curvature perturbation. Therefore, if these 
mixed derivatives are nonzero, the curvature perturbation is not conserved 
on large scales, even if the motion in field space is geodesic. 

When dealing with scalar fields, it is easier to work with the comoving 
curvature perturbation TZ, rather than with the curvature perturbation on 
uniform density hypersurfaces (. One can nonetheless write down the evo- 
lution equation for ( in the form (see [IQ] and ^Ti |12] for the exact and 
covariant form of this relation) 

^ = _i/^£i^_S, (60) 

p + p ^ ' 

where 5Pnad is the non-adiabatic pressure perturbation, defined as 

<5Pnad = 5P--5p, (61) 
P 
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and S is the scalar shear along comoving worldlines, which is given explicitly 
by 



H 



3H 



E - (B/a) + 



5q 



rC 



2p 



9{p + P) a?m 



(62) 



In our case, the non-adiabatic pressure perturbation can be expressed as 

^Prn 



SP 



nad 



QHXP 



X 



'l + ci)Pj-2XciP,xi 



;i + <)Pi 



2XctP 



XI 



(63) 

On sufficiently large scales, one can neglect the ffist term, proportional to 
the comoving energy density, and only the second term on the right hand 
side contributes to the non-adiabatic pressure perturbation. 



5 Two-field case 

We now specialize our formalism to the case where only two scalar fields are 
present. In this context, the entropy subspace is one-dimensional and the 
basis {e(, Cg} is completely specified. In order to make a direct comparison 
with the previous literature, we will use a more traditional notation, and 
replace the subscripts 1 and 2 by respectively a and s, so that 



5.1 Background equations 

The background equations of motion for the scalar fields can be decomposed 
into adiabatic and entropic equations. Defining 

a = V2X, (65) 

the adiabatic equation of motion can be written as 

a+(^3H + ^^a-^P^, = 0, P,. = eiPj. (66) 

Moreover, by using the decomposition 

Px = P,xxX + = P,xx ob^ P^xa a, (67) 
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the adiabatic equation of motion can also be rewritten as 



a 



P 



X 



P 



-3Ha 



.X 



{61 



The entropy part of the equations of motion gives the rate of change of the 
adiabatic basis vector in terms of the entropy basis vector ef: 



P,x^ 



-el. 



(69) 



This is the generahzation of the equation giving the time derivative of the 
angle 6 between the initial field basis and the adiabatic/entropy basis. This 
also implies that the non vanishing components of the matrix Z^n are given 
here by 

P. 



Px^ 



(70) 



To manipulate the equations of motion for the perturbations, the follow- 
ing identities will also be useful: 



P.r 



P. 



Pi 



aP 



X 



P P 



+ P^aaO- + P,XaCra, 



&P 



(71) 
(72) 



X 



5.2 Equations of motion for the perturbations 



Specializing our action ( l40l) to the present two-field case, one can express the 
action in terms of the quantities Q^- and Qs and easily derive their equations 
of motion. The adiabatic equation of motion can be written in the compact 
form 



Qa 



with 



and 



^^a 



3H 



P 



X 



p 



X 



\Hciy 



Qa + 
P,a 

aPx 



&p 



X 



cl)P, 



clP,XsCT^ 



&/H 



3H 



{a/Hy \ ja/Hy 
&/H J a/H 



(73) 
(74) 

(75) 
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The equation of motion for the entropy part is given by 



Qs 




^ En. (76) 



with 



PI 1 P P P 

'^^ = -i^ + r ^-2^n+2-p^, (77) 



where R denotes the Riemann scalar curvature of the field space. 

Using the relation (l58ll . this equation of motion can also be rewritten as 



Q, + (sH + ^) 4 + g + ,J + I) Q. = - - 4* . (78) 

which can be useful on large scales, when the right hand side can be neglected. 
In this limit, the above equation shows that the entropy perturbation Qg 
evolves independently of the adiabatic mode. 

It is clear from the previous equations that H quantifies the coupling be- 
tween the entropy and adiabatic modes. One also notices that the equations 
of motion are characterized by a friction term that differs from the usual 
Hubble friction term. The additional friction term for the entropy pertur- 
bation is the same as the one which appears in the background equation of 
motion, but different from the friction term for the adiabatic perturbation. 
Finally, as already emphasized for fields, we observe that the adiabatic 
mode propagates with the speed of sound Cg, while the isocurvature mode is 
characterized by a propagation with the usual speed of light. 

For a Lagrangian of the form 

P = X-Vi^'), (79) 

the mixing parameter reduces to H = —2Vs/<^ = —'^Zcrs- This shows that, in 
this particular case, the coupling between the adiabatic and entropy modes is 
directly related to the 'rotation' of the adiabatic/entropy basis in field space, 
i.e. to the bending of the background trajectory in field space. It is worth 
noting that this direct link no longer holds in the general case because of the 
term proportional to P^xs in S. When the field metric is fiat, Gjj = 6ij, one 
can introduce, as in [37] , the rotation angle bewteen the initial basis and the 
adiabatic/entropy basis. One thus finds that H = 26. For a Lagrangian of 
the form (17^ but with the field metric 

Gij d(t)^ dcj)-^ = dcj)'' + e^''^^) dx^ , (80) 
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one finds that E = 26 + b'a sin 6, where we see now that the additional 
term simply comes from the non-trivial covariant derivative in f lH^ . It is also 
easy to check that our equations ( jTSirfSll reduce to the results of [STJ HH] in 
these particular cases. Note that, in the first case, non-linear extensions of 
these equations have been obtained in [S], based on a covariant formalism 
similar to that of [ITl H^ . 

When one can neglect spatial gradients, the expression reduces to 




which implies that there exists a first integral for Qcr and that the second- 
order equation of motion (!73l) is not necessary in this limit. Indeed one can 
check that the large-scale limit of (!7S|) is a consequence of (IHTI) . 

To conclude this section, let us give the approximate form of the evolution 
equations for the perturbations when the spatial gradients can be neglected: 

(82) 




-^\Qs^O. (83) 



6 Quantum fluctuations 

We now consider the quantization of the system discussed in the previous 
section. In order to do so, we follow the usual procedure for single field infla- 
tion. The first step consists in introducing a new variable which is canonically 
normalized, with conformal time, and whose effective mass is time dependent 
because of the expansion of the Universe. In the present case, we have two 
degrees of freedom. By introducing the new fields 



aJP 



X 



Qa, Vs = a JPx Qs , (84) 

respectively for the adiabatic and entropy degrees of freedom, one can rewrite 
the second-order action in the form 



^(2) = \ I drd'k 



with 

a 



Cs ctPyc 



/2 I /2 n(- ' ; 2 2 2 ; 2 2 

h^laaVl + QssVs + 2QsaVaVs (85) 

[il + cl)P,s-cla'P,xs], (86) 
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and 



z" „ z' „ a" 



fio-cr — , ^s(j — — ^ , ^ss — ^ (87) 

z z a 



where we have introduced the two background-dependent functions 



z = \ , a = a^Px , (88) 

and used the conformal time r defined by r = / dt/a(t). 

For single field infiation with a standard kinetic term, z reduces to the 
usual function aa/H. In the general case considered here, the function z is 
still the ratio between and the comoving curvature perturbation defined 
in the previous section, i.e. 

v^ = zn. (89) 

The equations of motion derived from the action ( l85l) can be written in 
the compact form 

= 0. (90) 

= 0. (91) 

As mentioned already, the above system shows clearly that the adiabatic 
degree of freedom is sensitive to the sound horizon, while the entropy 
degree of freedom Vg is sensitive to the usual Hubble radius. In the general 
case where ^ is non-vanishing, the above system is coupled and we leave 
for a future investigation the analysis of the perturbations generated by this 
coupled system. In the following, we will just consider the simple case where 
^ can be neglected while the perturbations cross the Hubble radius and sound 
horizon. This is similar to the situation considered in [21] for two-field DBI 
infiation. 

When ^ = 0, the system is completely decoupled and one can analyse 
separately the adiabatic and entropy degrees of freedom. For completeness, 
let us repeat the standard analysis (see [26], or [39]). Each variable v, either 
Va- or Vg, becomes a quantum field which is decomposed as 

«(n ^) - / cl'* U^^^^y'-^ + at.Kr)e--^-n . (92) 




where the and a are creation and annihilation operators (there are two 
sets of creation and annihilitation operators, one for the adiabatic degree of 
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freedom, the other for the entropy one), which satisfy the usual commutation 
rules 



S{k - k') 



4' 4' 



0. 



(93) 



The action with = implies that the conjugate momenta for tv and Vg are 
respectively v'^ and v'^. Therefore, in either case, the canonical quantization 
for V and its conjugate momentum leads to the condition 



(94) 



which must be satisifed by the complex function f/c(r). 

For the adiabatic degree of freedom, the function tVfc(T) satisfies the 
equation of motion ( !90l) with ^ = 0, i.e. 



0. 



(95) 



In the slow-roll limit, and when Cg varies sufficiently slowly while the scale 
of interest crosses out the sound horizon, one can take the approximation 
z" /z ~ 2/r2, so that the general solution is known analytically within this 
approximation. Finally, we require that the solution on small scales behaves 
like the Minkowski vacuum. This leads to the solution 



1 



Vak 



-ikcs 



ylkcs 



kc.T 



(96) 



where the normalization is imposed by the condition ( 1941) . This implies that 
the power spectrum of the adiabatic fluctuations is given by 



k^ 
2^ 



Vak 



a^Px ^-K^CsPx' 



(97) 



where the quantities on the right hand side are evaluated at the sound horizon 
crossing. This can be translated into the power spectrum of the curvature 
perturbation IZ, 



n 



k^ \v<jk\'^ 
27r2 ^2 



{91 



where e = —H/H^. In the decoupled case, the adiabatic sector is equivalent 
to the single fleld A;-inflation scenario and the above expression coincides with 
the result given in [12]. A more reflned treatment, taking into account the 
next-to-leading-order corrections, can be found in the appendix of 
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Let us now consider the entropy degree of freedom. Wlien ^ = 0, tlie 
complex function ^^^(t) satisfies the equation of motion 

vl + (^^^ - ^ + = (99) 

In the slow-roll limit, e = —H/H^ is a small coefficient and its time derivative 
is at least second order in slow-roll. Since 

e = 2He{e + ^^^^\, (100) 



2HXP 



X , 



one must consider P x / {HPx) as first order in slow-roll. One can then write, 
around the time of Hubble-crossing {k = aH), 



v': + ev,-^U-^)v, = (101) 



with 



- 4 iJ2 + 2HP^x PxH^ 2H^ ^ ^ cl' am^P% ' ^ ' 



where we have neglected the square and the time derivative of P^x / {HP^x)- 
The solution of fllOip with the appropriate asymptotic behaviour is 

where H^'^ is the Hankel function of the first kind of order u. Note that we 
did not assume that the last three terms were small. If they are big, the 
effective mass is important and the entropy fluctuations are suppressed. The 
last term in particular shows that the effect of the bending of the background 
trajectory on the entropy mode is enhanced in the small sound speed limit. 
If the last three terms are small, i.e. in the limit Vg 3/2, one finds 

v,,^^e-''^(l-^]. (104) 



v/2A? V kr 

In this case, the entropy power spectrum, at Hubble crossing, is given by 

\Usk\ n. /incN 

T^Qs = o 9 2P - A 2P ' 

27r^ a^-T^x 4:'^' P,x 

where the quantities are evaluated at Hubble crossing. One sees that the 
ratio between the entropy and adiabatic power spectra, evaluated at their 
respective 'horizon' crossings, is therefore the speed of sound Cg. 
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7 Conclusion 



In this work, we have studied the hnear perturbations for a very large class of 
multi-field inflationary models, which includes most of the multi-field models 
studied so far. Our formalism should also be useful for studying the per- 
turbations in multi-field inflationary models that will be constructed in the 
future. In order to obtain the equations of motion for the perturbations, as 
well as their quantum initial conditions, we have calculated directly the ac- 
tion at second order in the linear perturbations. To facilitate the analysis of 
the equations of motion, we have decomposed the linear perturbations into 
one adiabatic mode, along the background inflationary trajectory, and — 1 
entropy modes. We have shown that a generic feature of these models is 
that the adiabatic and entropy modes obey coupled wave equations, which 
propagate with the speed of light for all entropy modes but with the speed 
of sound Cs for the adiabatic mode. 

We have also used this decomposition into adiabatic and entropy modes, 
to show that the (comoving) curvature perturbation, proportional to the adi- 
abatic perturbation, is sourced by a very specific combination of the entropy 
perturbations. In this combination, one finds a term which is non-zero when 
the background trajectory in field space is bent, or, more precisely, when it 
is non-geodesic (with respect to the field space metric G/j) but we have also 
found another term, proportional to P,xsi which is not present for models 
with canonical kinetic terms. 

Interestingly, despite (or maybe thanks to) the fact that we have consid- 
ered a much larger class of models than previous works, we have been able, in 
the two-field case, to write down the equations of motion in a very compact 
and simple form. 

Finally, it would be interesting to extend this work in several directions. 
The analysis of the quantum fluctuations when the adiabatic and entropy 
modes are coupled would be particularly interesting. It is more complicated 
than in the standard case, because the quantum-classical transition occurs 
in principle at different times for adiabatic and entropy modes with the same 
wave number k: at the usual Hubble crossing {k = aH) for the isocurvature 
modes, but at the crossing of the sound horizon [kcg = aH) for the adiabatic 
mode. 

In the present work, we have considered only linear perturbations. It 
would be interesting to go beyond the linear order and to study the non- 
Gaussianities that could be generated in this class of models. One can expect 
that the non-Gaussianities in this class of models would combine the types 
of non-Gaussianities exhibited in single-field DBI inflation and in multi-field 
inflation. 
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Appendix: multi-field DBI infiation 

An important example of infiation with non canonical kinetic term is the 
so-called DBI infiation model, where the infiaton field corresponds to the po- 
sition of a probe D3-brane moving in a warped background and its dynamics 
is described by a DBI action. Very recently, a few papers have considered 
a multifield extension of DBI infiation. Taking into account the angular po- 
sition of the brane in addition to its radial position, they started from an 
effective Lagrangian of the form 

P{X, 00 = -J (Vl - 2/A - l) - X = -^GjjV.c^'V^cP' (106) 

where / and V are functions of the scalar fields. The above Lagrangian be- 
longs to the class of models studied in the present work, and the results given 
in the main text can thus be applied to this particular case and compared 
with the results of [201 ED- 

The speed of sound derived from (11061) is given by 

c, = - 2/A = (107) 

where 7 is the analogue of the relativistic Lorentz factor. By taking the 
derivative of fll06p with respect to the kinetic term X, one sees that the DBI 
action is a very special case, which satisfies the identity 

Px = - = 1- (108) 

Since both Px and Cg appear very often in our equations for the perturba- 
tions, this special relation will bring simplifications. For example, the back- 
ground coefficients multiplying Q„ and Qg become, respectively, {3H + 3'y/'y) 
and {3H + 'y /'-/), in agreement with the equations given in [20] • Similarly, the 
coefficients of Ql and Q"^ in the action, which are respectively P^x / i^cj.) 
and a^P^xf^ reduce to a^/(2c^) and a^/(2cs) in the DBI obtained in 

We have also checked that, when specialized to the Lagrangian fll06p 
and the field metric (ISUjl . our equations for the perturbations agree with the 
equations given in [20]. Note, however, that the background coefficients in 
the equations of [20] are expressed in terms of the angle 6 between the initial 
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basis and the adiabatic/entropy basi^ in a form which makes them look 
singular when 6 = 0. In our case, we obtain directly the equations of motion 
with well-behaved coefficients. 

As far as the primordial power spectra are concerned, it is easy to verify 
that our expressions ( 197|) and (11051) reduce to the results of [21], when using 

(unzD- 

Let us finally note that the non standard contribution to the mixing 
coefficient, proportional to P^xs is non zero in general in this model since 

This term vanishes only when f^s vanishes. In the usual case where / is a 
function of the scalar field corresponding to the radial direction, / vanishes 
only if the background trajectory is strictly radial. 
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